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1 Introduction 



Non-relativistic QED (NRQED) |]J is becoming increasingly popular for QED bound 
state calculations |], ||, |3|, |J. It has the advantage over traditional Bethe-Salpeter equa- 
tions P| that the non-relativistic nature of the QED bound states is explicit, whereas 
relativistic and radiative corrections can be systematically incorporated by taking into 
account higher orders in the 1/m expansion and by calculating the matching coefficients 
at higher order in a respectively. However, the NRQED lagrangian still contains two 
dynamical scales, namely the typical relative momentum in the bound state p ~ ma 
and the bound state energy E ~ ma 2 , which implies that the terms in the lagrangian 
do not have a unique size. The leading size of each term is given by the next relevant 
scale p (soft) (except for the time derivative) and rules have been provided to estimate 
the subleading contributions due to the scale E (ultrasoft) 0. Nevertheless, it would be 
helpful for bound state calculations to have and effective field theory (EFT) where each 
term in the lagrangian had a well defined size. This EFT has proven to be quite elusive 
for some time fT]? H? H- 



In ref. |I(J we proposed potential NRQED (pNRQED) as such EFT and presented 
the form of its lagrangian for positronium. In ref. JTlJ we worked out pNRQED for 
Hydrogen-like atoms and reproduced the Lamb shift in a very straightforward way. It is 
the aim of this paper to discuss pNRQED for positronium in greater detail and to show 
that it also allows to reproduce the spectrum at order ma 5 , where all regions of momenta 
(hard, soft and ultrasoft) contribute, very efficiently. We also illustrate how dimensional 
regularization helps in that. 

pNRQED describes fermion-antifermion pairs with relative momentum of order p and 
energy of order E, and ultrasoft photons with energy and momentum of order E. This 
should be compared with NRQCD, which describes degrees of freedom (fermions and 
photons) with energy and momentum less than a certain cut-off fi such that E, p << jj << 
m. Formally speaking pNRQED has two UV cut-offs A x and A 2 , where E << A% << p 
is the cut-off for the energy of the fermions and for the energy and momentum of the 
ultrasoft photons, whereas p << A 2 << m is the cut-off for the relative momentum of 
the fermion-antifermion system. In principle, we have some freedom to choose the relative 
importance between Ai and A 2 . We choose A 2 /m << Ai, which guarantees that the UV 
behavior of the fermion propagators in pNRQED is that of the static ones. 
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pNRQED is obtained from NRQED by integrating out fermions and photons of ener- 
gies and momenta of order p and photons of energies of order E and momenta p[]. The 
pNRQED lagrangian obtained is local in time but non-local in space (i.e. it has potential 
terms) and contains ultrasoft photons only. The size of each term becomes explicit once 
the lagrangian is projected onto the one electron one positron subspace of the Fock space 
and written in terms of a wave function field. This is due to the fact that in the latter 
representation the ultrasoft photon fields can be multipole expanded about the center of 
mass. Moreover, the calculations in pNRQED are very close to those in nonrelativistic 
Quantum Mechanics. 

The practical way in which we integrate out degrees of freedom is by a matching 
procedure. We impose two fermion Green functions and four fermion Green functions 
(with an arbitrary number of ultrasoft photon legs) in NRQED be equal to those in 
pNRQED once both are expanded about the external fermion energies and (ultrasoft) 
photon energies and momenta. Dimensional regularization is used for both UV and IR 
divergences. Furthermore, we use static propagators, and hence the matching can be 
done to a given order in 1/m and a. This is justified because the fermion energies we are 
integrating out in loops are of order ma whereas the typical kinetic energy is 0(ma 2 ). 

The procedure above is similar to the matching between QED and NRQED as carried 
out in refs. |T^, [HJ. In that case scales ~ m are integrated out and the matching 



reduces to calculations in QED, where m is the only scale in the integrals. Here, the 
scales which are integrated out are ~ p and the matching reduces to calculations in 
NRQED where p is the only scale in the integrals. Hence, the potential terms in pNRQED 
play a role analogous to the Wilson coefficients in NRQED. Indeed, the former encode 
contributions due to physics at the scale p much in the same way as the latter do of 
physics at the scale m. At each matching step the non analytic behavior in the scale 
which is integrated out becomes explicit. 

We distribute the paper as follows. In sec. 2 we describe the matching procedure 
between NRQED and pNRQED. In sec. 3 we present the bound state calculation. Sec. 4 
is devoted to the conclusions and future prospects. In Appendix A gauge independence 



^ome authors p2| |lg|| like to distinguish potential photons, i.e. photons with k° ~ E and k ~ p from 
soft photons, i.e. photons with fc° ~ p and k ~ p. This distinction is quite irrelevant in our formulation 
since both potential and soft photons are integrated out at the same time when matching NRQED to 
pNRQED. 



4 



is checked at order ma A by calculating the matching to this order in the Feynman gauge. 
In Appendix B the pNRQED lagrangian for the unequal mass case is displayed. 



2 Matching NRQED to pNRQED 

The pieces of the NRQED lagrangian which are relevant to the calculation of the 
bound state energy at 0(ma 5 ) read 

, + f n D 2 D 4 a B [V-El 

Cnrqed = V [ <iD +- h--r + c F e- hc D e 



2m 8m 3 2m ' 8m 
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<r (DxE-ExD, , 

+ic S e - — - }ip + Xc, e -> -e 2.1 

Sm z 

m A m z 4 m z 

where ip is the Pauli spinor field that annihilates the fermion and XcQ is the Pauli spinor 
field that annihilates the antifermion. iD° = ido — eA°, iD = + eA on the fermion 
field. The bilinear lagrangian for \ c is equal to the ip one with the change e — > —e. 

This can be seen as follows. We draw all possible diagrams of the two fermion irre- 
ducible four fermion Green function, which can not be disconnected by cutting a photon 
line, such that s + r < 4, where s is the number of l/m factors in the diagram and r the 
number of explicit a H [K| [□]]. This rule can be easily justified if we take into account 



that the next relevant scale is p ~ ma and hence all l/m must be compensated for by ma 
until we reach dimensions of energy. For diagrams which can be disconnected by cutting 
a photon line the same rule applies but there is an extra suppression if n time deriva- 
tives act on this photon line. This is due to the fact that these time derivatives are only 
sensitive to the typical energy. The extra suppression factor is a n . Recall also that any 
diagram may have subleading contributions which appear from the analytic expansion of 
the external energy about zero. In the Coulomb gauge the diagrams fulfilling the above 
criteria are displayed in fig. 1 and fig. 2. Then both tree level and one loop diagrams are 
required. Notice also that cf, cd and cs are needed at order a and d s and d v at order a 2 . 
All these Wilson coefficients are gauge independent but depend on the renormalization 



2 We use here Xc — C\*, where C is the charge conjugation matrix, instead of \ as i n re f s - @; 
because it is more usual in nonrelativistic systems. 
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scheme for the UV divergences of NRQED. We shall use dimensional regularization with 



MS scheme. In this scheme the Wilson coefficients of the bilinear terms were given in 

15fl. They read 



14 and the ones for the four fermion terms in \M 




cf = 

CD = 

cs = 
d 2 = 

0U7T 

,/r 23 , tt\1 . . 

4 = ^ 1 -^i lo V + T- log2 + i 2j|' < 2 ' 2 > 

TTCt ( 2a (22 7T\] 

Hence our starting point is the lagrangian ( |2.1| ) with the Wilson coefficients ( |2.2j ). We 
wish to integrate out fermions and photons of energy and momenta ~ p and photons of 
energy ~ E and momentum ~ p. Then the effective theory we want to reach, namely 
pNRQED, contains only ultrasoft photons (of energy and momentum ~ E) and fermions 
of energy ~ E and momentum ~ p or less. Since we integrate out photon momenta ~ p 
but keep fermion momenta of this order, pNRQED contains terms non-local in space, 
namely potential terms. This is not a problem for a non-relativistic EFT. 

The practical way to obtain pNRQED is by enforcing two and four fermion Green 
functions with arbitrary ultrasoft external photons to be equal to those of NRQED once 
we expand about zero the energy in the external electron legs and the energy and momenta 
of the ultrasoft photon legs. This may produce IR divergences which are regulated in 
dimensional regularization, in the same way as the UV divergences are. Since the IR 
behavior of NRQED and pNRQED is the same, these IR divergences will cancel out in 
the matching. The UV divergences of NRQED must be renormalized in the MS if we 
want to use the matching coefficients (|2.2|). We still have a choice in the renormalization 



scheme of pNRQED. However, it is most advantageous to use again MS. Indeed, in this 



scheme we can blindly use MS for any divergence regardless it is UV or IR in the matching 
calculation. For the UV divergences of NRQED and pNRQED it is just the scheme we 
choose, and for the IR divergences it is irrelevant as long as we use the same treatment 
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in both theories, since the IR behavior is the same. This allows to put integrals with no 
scale equal to zero and, as we will see later on, reduces the matching to a calculation in 
NRQED only. 

Notice that we demand off-shell Green functions in NRQED and pNRQED to be 
equal and not on-shell Green functions as it is costume in the matching from QED to 
NRQED Jl], pL4|| . This is due to the fact that we are eventually interested in bound 
states, and particles in a bound state are typically off-shell. The equations of motion of 
pNRQED (with potential terms included), or local field redefinitions, may be consistently 
used later on to remove time derivatives in higher order terms and write the pNRQED 
lagrangian in a standard form, within the philosophy advocated in ref. |16| (see also 
|17|). We have checked in Appendix A that this procedure produces gauge independent 
results at O(ma 4 )0. In fact the same argument applies to the matching between QED and 
NRQED, which accordingly should also be carried out off-shell. However, in that case, 
at lower orders there is no difference between doing the matching on-shell and replacing 
derivatives by covariant derivatives, in order to enforce gauge invariance, and doing the 
matching off-shell and consistently using the equations of motion or local field redefinitions 
to get rid of time derivative terms [14], [16| . 

The remaining important ingredient to carry out the matching efficiently is the use 
of static (HQET) propagators for the fermions. This has been completely justified in the 
matching between QED and NRQED [14], [15[] , since both energy and momentum in this 
theory cannot exceed the same cut-off A which is smaller than m. Hence the UV behavior 
of the fermion propagator in NRQED is always dominated by the energy. This fact is not 
automatically implemented in dimensional regularization. When dimensional regulariza- 
tion is used, the correct UV behavior of NRQED is only obtained when expanding about 
the static propagator. 

In pNRQED we have a certain choice for the UV cut-offs for the fermion energy and 
momentum. We shall choose Ai << p for the energy and A2 << m for the momentum in 
such a way that K\jm << A x . The proper way to implement this condition in dimensional 

3 However, there is still some freedom in the choice of the wave function field, to be introduced later 
on, due to time independent unitary transformations which commute with the leading terms in the 
pNRQED lagrangian. Therefore, in general, it is not to be expected that the standard form of the 
pNRQED lagrangian calculated with different gauges coincide, but only to be related by one such unitary 
transformation. This explains, for instance, why the potential presented in [is[| is different from ours but 
leads to the same physics. 
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regularization is again by expanding the fermion propagator in pNRQED about the static 
propagator. 

Now we are in a position to prove that no pNRQED diagram containing a loop con- 
tributes to the matching calculation. 

Consider first the two fermion Green function with an arbitrary number of ultrasoft 
legs. For potential terms to contribute we need at least a four fermion Green function and 
hence we only have to care about ultrasoft photons. If we input a momentum ~ p in the 
fermion line this momentum cannot flow out through any external ultrasoft photon line 
(by definition of ultrasoft). Then it must flow through the fermion line, which is a series 
of static propagators insensible to the momentum flowing. Hence upon expanding about 
external fermion energies and external energies and momenta of the ultrasoft photons there 
is no scale in any integral and therefore any loop contribution vanishes. In fact, exactly 
the same argumentation can be used for the NRQED calculation. Then we conclude that 
the terms bilinear in fermions are exactly the same in NRQED and pNRQED. However, 
we have to keep in mind that the latter (by definition) must be understood as containing 
ultrasoft photons only. 

Consider next the four fermion Green function in pNRQED containing several po- 
tential terms but no US photon. Since no energy can flow through the potentials and 
the static propagators are insensitive to the momentum, upon expanding about the US 
external energy, the integrals over internal energies have no scale. However, these inte- 
grals have IR (pinch) singularities which are not regulated by standard DR. We shall take 
the additional prescription of putting them to zero. Since the IR behavior of pNRQED 
and NRQED is the same, the same kind of integrals appear in the NRQED calculation. 
If we put them consistently to zero we obtain the correct potential terms, which play 
a role similar to the Wilson coefficients in the matching between QED and NRQED. It 
is important to keep in mind that the Wilson coefficients compensate for the different 
UV behavior of the effective theory (pNRQED) with respect to that of the 'fundamental' 
theory (NRQED). Hence they are not sensitive to the details of the IR behavior, which 
legitimates the prescription above. Then any loop diagram in pNRQED with no US pho- 
tons can be put to zero. This still holds if an arbitrary number of US photon lines are 
included in the diagram. Indeed, any potential line in the diagram now may also con- 
tain US momenta from the photon lines. These, however, can be expanded about zero 
since they are (by definition) much smaller than the momentum transfer in the potential. 
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Hence the integrals over US photon energies and momenta contain no scale (again upon 
expanding the US external energy in the fermion static propagators) and can also be put 
to zero. 

In summary, we can directly identify the potential terms from a calculation in NRQED. 
We would like to stress again the similarity in the procedure with the matching between 
QED and NRQED as carried out in refs. |L4], |l^]. The potential terms in pNRQED play 
the role of Wilson coefficients in the matching procedure. 

The four fermion terms appearing in the pNRQED lagrangian typically have the form 4 
L pot = -J d 3 ^d 3 ^(t, xi)xJ(t, x 2 )U(x, p, (Ti, (T 2 )Xc{t, x 2 )V(t, xi) (2.3) 

where x = x x — x 2 , p = — «V X and Si = 0"i/2, s 2 = <x 2 /2 act on the fermion and 
antifermion, respectively (the spin fermion and antifermion indices are contracted with 
the potential indices, which are not explicitly displayed). K(x, p, cr 1; cr 2 ) may also be 
written as an expansion of the type 

V = V {0) + U (1) + U (2) + y (3) + ■ ■ • , (2.4) 

where (V^) ~ ma n . Our results are exact for the four first terms of this expansion. 

We obtain from the tree level diagrams of fig. 1 (V represents the Fourier transform 
of V) 
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In principle (2.3) could also depend on the total momentum P = — zVx, with X 



ultrasoft photons, but these effects can be neglected to the accuracy we are working at. 
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v£ = ^ U • * - . (2-11) 



??? 



k 2 



v£?e = (d a + 34) - 2^S 2 , (2. 12) 

where S = sj + S2, and for the one loop diagrams of fig. 2. 

^i=^flog^-?log2 + ll , (2-13) 



^- = S( 1 °V + 2l0g2 - 1 ) • (2 - M) 

The \i dependence of fig. 2b, c is of IR origin and will eventually cancel with US contri- 
butions. However, the /i dependence of fig. 2a is of UV origin and cancels exactly with 
the ji dependence of d s . Recall that there is an additional // dependence in cd which will 
also cancel against US contributions. Upon Fourier transforming and putting together 
the above results we obtain 

V = V {0) +5V J (2.15) 



where 



= ~ (2.16) 



and 



a 1 / 2 1 



j(3) ( x ) ( ( o2\ j 1C j " 2 7a \ A 7a 2 1 

H 5— hra c D - 2Cp } + d s + 3d v - 167rad 2 + — — log/i - - -reg-— 7 

m A V 06) birm z \x\ 6 



, -S 2 (vra-4 - 24) + t^^L • S (2c s + Ac F ) + ^f — S 12 (x) , 
rrr \ 3 / 4m z xr 4m / xr 



where Si 2 (x) = (— <Ti • cr 2 + 3<T! • xcr 2 • x) and (see |18| for more details on the Fourier 
transform) 



The pNRQED lagrangian now reads 



- J d 3 X!CpX.2^(t, *l)xl(t, X 2 )V(X, p, <7i, (T 2 )Xc(t, X 2 )^(t, X X ) , (2.19) 
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where the photons are ultrasoft. 

In order to make explicit the size of each term in (|2.19 ) it is convenient to project 
pNRQED to the one electron one positron subspace (this can be easily done at the hamil- 
tonian level). This subspace is spanned by 

J d 3 x l0 ? 3 x 2 ^( Xl , x 2 )^ t (x 1 )xl (x 2 ) |0) , (2.20) 

where |0) is the subspace of the Fock space containing zero electrons and positrons but 
an arbitrary number of ultrasoft photons. 

Then the dynamics of the wave function field is described by the lagrangian 

f ( V 2 V 2 V 4 

L p nrqed = / t^xi^xz^x!, x 2 , t) [ id + — + — ^ + —± - ex l diA (t, X) 
J \ m Am Am 6 

-2ie A(t ' X) ' Vx - V(x, p, <t u <r 2 ) Wi, x 2 , t) - f d?*-F^\ (2.21) 
m J J 4 

where we have made precise that the remaining photon fields are ultrasoft by multipole 
expanding them about the center of massQ Furthermore, gauge invariance at any order 
in the multipole expansion can be made explicit by introducing 

<p(xi, x 2 , t) = P[e ie £* A dx ]5(x, X, t) . (2.22) 

Then, the gauge transformations of the above wave function fields are 

<p(xi, x 2 , t) -> </(xi, t)^(xi, x 2 , t)5f _1 (x 2 , t) , (2.23) 

S(x,X,t) -> S(x,X,t). 

We finally obtain 

LpNRQED = I dVXdtS\x,X,t){id - + - (2.24) 
J [ m 4m ci 4m 

-V(x, p, tr 1; <r 2 ) + ex • E(X, t) js(x, X, i) - | d 3 xi^F^ , 

which is explicitly gauge invariant. Moreover the size of each term is unique and can be 
evaluated as follows. Each relative momentum <9 X , and inverse relative coordinate |x| _1 
have a size ~ ma. Each US photon field, derivatives acting on it, the time derivative and 



3 As expected for a chargeless particle the covariant derivatives for P are the ordinary ones. 
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the center of mass momentum <9x (in the rest frame, when entering in recoil corrections 
due to the virtual emission of US photons) on the wave function field have a size ~ ma 2 . 
Formula (2.23) has already been presented in |L0| (except for a numerical factor in the 
potential which was left arbitrary). We shall use it to calculate the spectrum at 0(ma 5 ) 
in the next section. 

The gauge independence of the matching calculation is checked at 0(ma A ) in Appendix 
A. The pNRQED Lagrangian for the unequal mass case can be built with no further 
difficulty The result is displayed in Appendix B. 



3 Bound state calculation in pNRQED 

In order to find the corrections to the bound state energy of a state with principal 
quantum number n we consider the following Green function (we will follow a procedure 
and notation similar to ref. ||1 1|| ) 

IT(g,x) := J dx°dXe igx ° < T{^(0)^ t (x, X, t)} > (3.1) 

when q —>■ E n , where E n is the energy of the leading hamiltonian 

V 2 a 

h = fr. (3.2) 

m |x| 

The integral over X fixes the center of mass momentum P to zero. We write 

U(q, x) = A " + 6An „ A - ± 5An + _J*2_ tf ^_L_ • (3.3) 
q - [E n + 6 E n ) q-E n q - E n q - E n 

The contribution to 5E n coming from the correction to the potential and the kinetic 
energy read 

ma 4 f 3 1 

IT \n 3 (2l + 1) + ^ 



S v E n = (nlj\5V\nlj) 



2a Sio ( m [in \1 ( " 1 n - 1 

-log— + 71og 61og2 + — -7 2^7 + ^ — 

Sir n 6 \ n ma 5 xt^i ^ ^ n 

7a 1- 5io 1 



3vr n 3 /(/ + 1)(2/ + 1) 



(3.4) 



US l0 S sl ( 3a f 32 \\ (1 - 6 l0 )5 sl ' 

+ T IT I 1 + {-g - 21og2j } + m + m + 1)w3 Cj,i 
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Am 3 



(nlj\V 4 \nlj) 



ma 4 3(1 + 1/2) -An 
~8 4n 4 (2/ + 1) 



(3.5) 



where 



3,1 



l + l 

21- 1 

a 

-2 

I * 
21 + 3 



a 



2(31 - 1) + -(Al - 1 



7T 



, 3 
■ J 



(3.6) 



a 



2(3l + A) + -(Al + 5) 



7T 



3=1 + 1. 



There is also a contribution from a virtual exchange of an ultrasoft photon corre- 
sponding to the diagram in fig. 3, which has already been evaluated in dimensional 
regularization for the Hydrogen-like atom |TTJ (here the calculation is identical but us- 
ing the reduced mass). Notice that the MS scheme has to be used in the calculation. 
Since (|2.24j ) is gauge invariant we can use any gauge to calculate this contribution. Still 
the Coulomb gauge continues to be advantageous, since in this gauge A can only con- 
tribute to tadpoles which can be safely put to zero in dimensional regularization. This 
contribution reads 



6 us E n = -§° £ \(n\ P |m)| 2 (E n - E m ) (log 
3tt ^ m \ 



E„ — E n 



log2 + - 



(3.7) 



ma 
3im 3 



5, 
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log 



+ (l-<J,o)log 



m 2 a 4 



where the last equation implicitly defines (E n j). Once we add all these contributions the 
final result reads 

.,4 

6E n /„• = 5 v E n + 5 K E n + 6 us E n 



ma 
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8 { n 3 (2l + l) 8n 4 

' 49 / n 1 

9 logo; + 7 log 7i + 8 log i2(n,Z) -14 log 2 7 V- + 
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3tt 
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log R(n,l) + 



3a f 32 



1 



16/(/ + l)(2/ + l) 



(31 



1 + — -— -2 log 2 



7tt 



where log R(n, I) = log 2 ^'2^ is called the Bethe logarithm. For the 0(ma 5 ) contribution 
we find agreement with ref. |1| for the spin independent piece and for the 5;o^si piece, 
while for the (1 — <5; )<5 s i piece we find agreement with ref. || (this last piece could also 
be obtained from results of ref. fll9|). 
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1(21 + !)(/ + l)n : 
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We can also easily obtain the full decay width at lowest order. It reads 

r B = ^*io(l-i,i)+ £ ^-a\(n\^-\m)\ 2 (E n -E m ). (3.9) 



4 Discussion and Conclusions 

We have seen that pNRQED correctly reproduces the positronium spectrum at 0(ma 5 ). 
This is a non-trivial check of the ideas behind this EFT since at this order all regions of 
momenta (hard, soft and ultrasoft) contribute to the energy. 

We would like to emphasize that the procedure we propose for higher order corrections 
to the positronium (and other QED bound states) is totally systematic. It uses two EFTs, 
namely NRQED and pNRQED. Both the matching from QED to NRQED and from 
NRQED to pNRQED can be done order by order in 1/m and a, and static propagators 
for the fermions can be used. This together with the use of dimensional regularization 
simplifies a lot the calculations. The actual bound state calculation is done at the level of 
pNRQED for the wave function field and is very similar to a standard quantum mechanical 
calculation. The only difference being that the wave function field couples to US photons 
in a field theoretical fashion. 

We believe that the clarity and simplicity of this formalism will allow to carry out 
higher order bound state calculations in QED very efficiently In order to illustrate this 
point let us pose ourselves the calculation of the positronium spectrum at 0(ma 6 ) and see 
the extra calculations required in order to obtain the pNRQED lagrangian at this order. 
Clearly all contributions that we obtain at 0(ma 5 ) which are multiplied by a Wilson 
coefficient, will give a contribution at 0(ma 6 ) by just calculating the Wilson coefficient 
to next order in a. This requires the matching from QED to NRQED at two loops. The 
relativistic correction to the kinetic energy 0(l/m 5 ) in NRQED should be kept. Terms 
0(l/m 4 ) in the NRQED lagrangian would now contribute to the potential but they do 
so only at tree level. Hence their Wilson coefficients are only necessary at tree level. 



Those for the bilinear terms may be obtained from []20j whereas those for the four fermion 
operators from pifl . Terms 0(l/m 3 ) also contribute at tree level and may contribute at 



one loop. In either case the Wilson coefficients are only needed at tree level which are 
known flij| . In addition 0(l/m 2 ) terms in the NRQED lagrangian would now contribute 
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to the potential at one loop and 0(l/m) terms at two loops. It is also easy to see by 
inspecting the next order terms of the multipole expansion in the pNRQED lagrangian 
that, due to angular momentum conservation, no contribution of US photons arises at 
0(ma 6 ). There would only be a new qualitative feature, namely that time derivatives 
multiplying potential terms would arise (for instance, from the expansion of the energy in 
the one transverse photon exchange at tree level). This time derivatives can be disposed 
of by using the equations of motion in pNRQED (now with potential terms included) as 
it has been done in Appendix A, according to the philosophy of ref. [[U| (see also []17|). 
This calculation would produce an independent check of the existing results obtained very 



recently in ref. 22 



It is important to be able to calculate systematically higher orders in QED as a test to 
the Standard Model in a sector where QCD does not play any relevant role. Any signal of 
new physics here should be much clearer than in the hadronic sector as, for instance, in 
the orthopositronium decay where there seem to be some difficulties to explain the data 
(see [E3| and references there in)F|. 



Finally, we would like to stress that the idea of separating the calculation of the 
binding energy (or any other observable) of a non-relativistic bound state system in three 
stages, namely (i) integrating out the hard scale, (ii) integrating out the soft scale and (iii) 
calculating the bound state energy when only ultrasoft degrees of freedom are present, is 
neither confined to positronium nor to the QED realm. Indeed, we have already shown 
in ref. |ll| that it can be applied to Hydrogen-like atoms, and it should be easy to work 
out pNRQED for muonium, di-muonium, and other two body QED bound states. In 
particular, pionium, a QED bound state which, however, decays strongly, has received 
considerable attention lately f25 |. Its decay width turns out to be proportional to the 
pion scattering length, which is an essential input to fix the parameters of the chiral 



lagrangian ||26|| . In order to extract the scattering length neatly from the experimental 
data a good control on the electromagnetic corrections is necessary. pNRQED for pionium 
can definitely help in that goal. Beyond QED, heavy quarkonium systems also form non- 
relativistic bound states. We have already proposed that potential NRQCD (pNRQCD), 
an EFT for NRQCD analogous to pNRQED for NRQED, should be useful to study 
these systems [|Kj. The techniques presented here may also help in the understanding of 
nucleon-nucleon bound states from the heavy baryon chiral lagrangian, which have also 

6 Only one experiment |24[] seems to be compatible with theoretical predictions. 
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received quite some attention during the last years p7 |. 
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Appendix A : 0(ma A ) matching in the Feynman gauge 

In this Appendix we check that the pNRQED lagrangian is gauge independent at 
0(ma 4 ) once it is written in the standard form. By the standard form we understand 
that all time derivatives in higher order terms are disposed of by the use of equations of 
motion, or alternatively by local field redefinitions [H, ITBl (see also |T7j). In order to do 



so the matching calculation is redone in the Feynman gauge. 

The main difference between the Coulomb gauge and the Feynman gauge, as far as 
the matching of the four point Green function is concerned, is that in the former loop 
diagrams involving A only can be set to zero in NRQED because there is no scale for the 
integration over the energy, whereas in the latter they must be kept because the poles of 
the Ao propagator now relate energy and momentum. 

The counting in section 2 implies now that the following extra diagrams must be 
considered: (i) the one loop diagrams of fig. 4 (0(ma 3 )), (ii) the same diagrams with a 
p 2 /2m insertion in either fermionic line (0(ma 4 )), (iii) the same diagrams with external 
energy insertions arising from the expansion of the propagators about zero external energy 
(0(ma 4 )), and (iv) two loop diagrams involving Aq propagators only (0(ma 4 )). 

It is easy to see that diagrams (i) cancel each other. Diagrams (ii) and (iii) vanish 
individually, essentially because they have an odd number of static propagators. Diagrams 
(iv) have already been seen to cancel in ref . 2l| . Then we are left with the same diagrams 



we had in the Coulomb gauge (all diagrams in fig. 1, except (d) which is 0(m« 5 )), but 
now they must be calculated in the Feynman gauge. 

In fact, all the diagrams give the same result except fig. la and fig. le. The latter 
now reads 

.2 



^(Feynman) = ^(Coulomb) - — 2 [P J 1 . (A.l) 
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Fig la now receives a contribution due to the expansion of the external energies about 
zero in the A propagator, which we depicted in (fig. 5). It reads 

^(Feynman) = -Ana (^) , (A.2) 

where k° = E x - E[. E l and E[ give rise to time derivatives in the pNRQED lagrangian. 

L"* (Feynman) = -J d 3 ^d 3 ^ 2 d 2 (^ty(f,xi)) / 7^8 e *" X ^xJxc(*,x a ) 

= / d 3 Xl rf 3 x 2 9 (vVCt.xO) / -^e*-*ld (xkc(t,x 2 )) (A.3) 

We can get rid of these time derivatives by using the equations of motion. Notice however 
that now potential terms enter in the equations of motion. Their explicit inclusion can 
be avoided in this case by using the continuity equation in the last equality of ( |A.3| ) 

!+V.J = 0, (A.4) 



where p = ip^ip and 



i 

2m 



ip^Vip - {Vip^ip] . (A.5) 



( |A.2| ) can now be written as 



CeT(Feynman) = +— V ^ ' , (A.6) 

which just cancels the extra contribution in (|A.1| ). We have then proved that the pNRQED 
lagrangian written in the standard form (i.e. with no time derivatives) at 0(ma 4 ) is 
exactly the same in the Coulomb and Feynman gauges. Notice that it has been crucial to 
write the time derivatives in a symmetric fashion in order to use the continuity equation. 
The naive use of the on-shell condition k° = p 2 /2m — p /2 /2m in (A.2) leads to incorrect 
results. 



Appendix B : pNRQED lagrangian for the unequal 

mass case 

Here we display the lagrangian relevant for the calculation of the mass to 0(ma 5 ) for 
the unequal mass case (we assume mi, m 2 >> p >> E). The charge of each particle has 
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opposite sign. 



= /^X^(x,X, t ){ ia „ - £ + j£ + - ^ (B.1) 

-V(x,p,a 1 , CT2 )+ex.E(X,i)Js(x,X,«)- | tftcii^f" , 

where M = mi + m 2 , /ii2 = m^+m, 2 ' x anc ^ an< ^ P anc ^ ^ are ^ e re l a ^i ve an d center of 
mass coordinate and momentum respectively. The potential now reads 

^ = - n - ~ 1 r P 2 + -x • (x • p p (B.2) 

|x| lvn\vn2 |x| \ x 2 / 

+ 5( 3 )(x) 



7a 2 



rra I gg^l±^i _ 24 ) + ds + 3^ - 16^ 2 f m ' + mi > | + ^ 



3 



log//' 



7a 2 1 5 (3) (x) c2 / 4 2 

reg — n H b ixa-c F 

67rmim2 Ixl mim 2 V 3 



acp 1 ^ acs 1 ^ /s 2 m 2 + s 1 m2 N \ acf* 1 

m\m,2 Ixl 3 2m\m2 Ixl 3 \ mim 2 / 4mim 2 |x| : 



where 



a / 4 , m. 2 ' 

71 \3 \l l 



and now (i s and d„ read (see 



tt 2 771? 

^ = — 9 9 -m 1 m 2 log — „- . (B.5) 

mf — m 2 ?n 2 

The lagrangian (B.2) must be corrected if there are charged particles of masses m^, 
i = 3,4.., similar or smaller than m := max{mi,m 2 }. Each particle of mass such 
that m > rrii » fiuoc gives an extra contribution 1/mf multiplying to g? 2 in (B.2). Each 
particle of mass [ivioi > >> //i 2 a 2 gives extra non-trivial contributions to the potential 
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Captionl. The nonzero relevant diagrams for the matching at tree level in the 
Coulomb gauge. The dashed and zigzag lines represent the A and A fields respec- 
tively, while the continuous lines represent the fermion and antifermion fields. The first 
diagram is the Coulomb potential. For the A the circle is the vertex proportional to c D , 
the square to cs (spin dependent) and the dashed dot to d 2 (the vacuum polarization), 
while for A the square is the vertex proportional to cp and the other vertex appear from 
the covariant derivative in the kinetic term. The last diagram is proportional to d s and 
d v . The symmetric diagrams are not displayed. 

Caption2. The nonzero relevant diagrams for the matching at one loop in the 
Coulomb gauge. The dashed and zigzag lines represent the A and A fields respectively, 
while the continuous lines represent the fermion and antifermion. The interactions for 
A are the ones which appear from the covariant space derivatives in the kinetic term, 
while for Aq comes from the covariant time derivative. The symmetric diagrams are not 
displayed. 

Caption3. The thick line and wavy line are the positronium and the transverse 
ultrasoft photon propagators respectively. 

Caption4. The nonzero relevant diagrams for the matching at one loop in the Feyn- 
man gauge setting E to zero and with no 1/m insertions. 

Caption5. Correction to the A propagator due to energy insertions in the Feynman 
gauge. 
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